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The combination of the unique properties of graphene with spin polarization and magnetism for 
the design of new spintronic concepts and devices has been hampered by the smaU Coulomb interac- 
tion and the tiny spin-orbit coupling of carbon in pristine graphene. Such device concepts would take 
advantage of the control of the spin degree of freedom utilizing the widely available electric fields 
in electronics or of topological transport mechanisms such as the conjectured quantum anomalous 
Hall effect. Here we show, using first-principles methods, that 5d transition-metal (TM) adatoms 
deposited on graphene display remarkable magnetic properties. All considered TM adatoms possess 
significant spin moments with colossal magnetocrystalline anisotropy energies as large as 50 meV per 
TM atom. We reveal that the magneto-electric response of deposited TM atoms is extremely strong 
and in some cases offers even the possibility to switch the spontaneous magnetization direction by 
a moderate external electric field. We predict that an electrically tunable quantum anomalous Hall 
effect can be observed in this type of hybrid materials. 



Spin-orbit interaction, which couples the spin degree 
of freedom of electrons to their orbital motion in the 
lattice, leads to many prominent physical phenomena, 
such as the antisymmetric exchange interaction 1 , the 
colossal magnetic anisotropy [5], or the anomalous Hall 
effect [3] in conventional ferromagnets. It is also the 
key interaction in the newly-found quantum topologi- 
cal phase in topological insulators, where the quantum 
spin Hall effect was observed experimentally [4], and the 
quantum anomalous Hall effect (QAHE) was predicted 
to exist [51 ini ■ Since the orbital motion can be manip- 
ulated with external electric fields, spin-orbit coupling 
(SOC) opens a route to dissipationless transport and to 
electrical control of magnetic properties [THS] , playing a 
crucial role in future spintronics applications. 

Owing to its strong spin-orbit coupling, heavy Ad and 
5d TMs display fascinating physical properties for desir- 
able spintronic applications, especially when combined 
with non-vanishing magnetization. However, magnetism 
of 5d TMs proved difficult to achieve due to their more 
delocalized valence d wavefunctions and the smaller intra- 
atomic exchange integrals compared to 3d TMs such as 
Fe, Co, or Ni. Low-dimensional structures, such as mono- 
layers [TO: , atomic wires [TT] , and clusters [T2] , were pro- 
posed to facilitate the stability of 5d magnetism by reduc- 
ing the d-d hybridization. Nevertheless, when deposited 
on substrates of noble or late transition metals, such as 
Cu, Ag, Au, and Pt, interdiffusion at interfaces is in- 
evitable, and strong hybridization between the d states 
of the substrate and of the adatoms is also destructive 
for 5d magnetism. From this point of view, using sp sub- 
strates is more promising. In fact, 4c? ferromagnetism 
was first observed in a Ru monolayer deposited on the 
graphite (0001) surface [13], which is close in its chemi- 
cal and physical properties to graphene. 



Since the day graphene was isolated and produced as 
a two-dimensional material I 1 . it abruptly altered the 
research direction of material science p[5] with the aim 
of exploring its fascinating transport properties. In a 
sense, graphene serves as a prototype of topological in- 
sulators [IF. For instance, the Berry phase of tt of elec- 
tronic states in graphene induces a half-integer quantum 
Hall effect [171 HH] > while the existence of the quantum 
spin Hall effect was first suggested for pure graphene 
when SOC is taken into account [19) . One common 
feature for those transport properties is the non-trivial 
topological origin, resulting in dissipationless charge or 
spin current carried by edge states with conductivity 
quantized in units of e'^/h. However, from the appli- 
cation point of view, a large external magnetic field is 
required to obtain the quantum Hall effect, and the spin- 
degeneracy in the quantum spin Hall effect makes it hard 
to manipulate the spin degree of freedom by controlling 
external fields. To avoid such constraints while keeping 
the benefit of topologically protected quantized topolog- 
ical transport, the long-sought QAHE is a perfect solu- 
tion. The essence of the QAHE lies in the quantization 
of the transverse charge conductivity in a material with 
intrinsic non-vanishing magnetization. The fact that the 
magnetization in ferromagnets can be much easier han- 
dled experimentally than large magnetic fields makes the 
QAHE extremely attractive for applications in spintron- 
ics and quantum information. However, at present the 
QAHE is merely a generic theoretical concept for mag- 
netically doped topological insulators [3 [6|. Recently, 
Qiao et al. suggested that the QAHE could also occur at 
comparatively low temperatures in graphene decorated 
with Fe adatoms [20] . 

Here, we demonstrate based on first-principles theory 
that 5d TMs deposited on graphene are strongly mag- 
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netic, provide colossal magneto-crystalline anisotropy en- 
ergies and exhibit topologically non-trivial band gaps due 
to very strong spin-orbit interaction. A generic represen- 
tative of this hybrid class of materials has a magneto- 
crystalline anisotropy energy of 10—30 meV per TM and 
a quantum anomalous Hall gap in its electronic spectrum 
with the size of 20—80 meV. In connection with a large 
magneto-electric response of the deposited adatoms, we 
predict that within this class of systems an electrically 
tunable QAHE at room temperature could be achieved 
experimentally. 

We have performed first-principles calculations of 5d 
TM (Hf, Ta, W, Re, Os, Ir, and Pt) adatoms on graphene 
in 2x2, 3x3, and 4x4 supercell geometries corresponding 
to the deposition density of 4.7, 2.1 and 1.2 atoms/nm^, 
respectively (see the Method part for details). Through- 
out this work the TM atoms are placed at the hollow 
sites of graphene. According to previous theoretical stud- 
ies [2T], it is the most favorable absorption site with the 
exceptions of Pt [52] and Ir |23j, for which the bridge 
site is more preferable. Among all systems studied we 
selected one prototype system, W on graphene, that we 
discuss in more detail. In the 4x4 geometry, the opti- 
mized distance between W adatoms on the hollow site 
and the C plane is about 1.74 A, and the hollow site is 
about 0.14 eV (0.41 eV) per TM atom lower in energy 
than the bridge (top) site, in agreement with previous 
observations [24]. Applying an electric field, the relaxed 
positions change by at most 0.01 Bohr radii, and the hol- 
low site remains the preferred adsorption site. 

Turning now to magnetism, and looking first at TM 
atoms adsorbed in the 4x4 supercell, we find that most 
of the 5d TM atoms (except for Pt) are magnetic with 
sizable magnetic moments (see Fig. la) ranging between 
0.5 and 2fiB- This is consistent with the large magnetiza- 
tion energy calculated, defined as the total energy differ- 
ence between the nonmagnetic and ferromagnetic state, 
reaching e.g. for Ta and W values as high as 0.56 eV 
and 0.34 eV, respectively. Surprisingly, the magnetic mo- 
ments across the series of the rather isolated TM atoms 
do not follow Hund's first rule, i.e. the magnetic mo- 
ments do not increase in steps of IfiB from Hf to Re, 
reaching a maximum at the center of the TM series and 
decreasing again from Re to Pt. On the contrary we 
find the minimal spin moment for Re. This demonstrates 
the role of the hybridization of the TM d-orbitals with 
graphene. Increasing the deposition density studying a 
3x3 supercell, the results remain basically unaltered, in- 
dicating that the d-d hybridization between neighboring 
adatoms is week and that the physics is determined by 
the local TM-graphene bonding. This changes in the 
case of the 2x2 geometry, when the adatoms are in close 
vicinity to each other and the magnetism survives only 
for Ta, W, Re, and Ir atoms. 

Unique to the 5c? TMs is the property that the SOC and 
the intra-atomic exchange are of the same magnitude. 



From this we anticipate a significant influence of SOC on 
the electronic properties, leading to strong anisotropies 
of the calculated quantities upon changing the direction 
of magnetization M relative to the graphene plane. 

The most impressive manifestation of SOC in the TM- 
graphene hybrid materials are the colossal values of the 
magneto-crystalline anisotropy energy (MAE), defined as 
the total-energy difference between the magnetic states 
with spin moments aligned in the plane and out of the 
graphene plane. The MAE presents one of the most fun- 
damental quantities of any magnetic system as its sign 
defines the easy magnetization axis and its magnitude 
gives an estimate on the stability of the magnetization 
with respect to temperature fiuctuations — a key issue 
for magnetic storage materials. A large MAE makes the 
magnetization very stable but also difficult to manipu- 
late. 

The values of the MAE obtained in 3x3 and 4x4 ge- 
ometry, presented in Fig. lb, lie in the range of 10 to 
50 meV per TM. These are orders of magnitude larger 
than the values for 3d TMs, irrespective of whether they 
are arranged in bulk, thin films or multilayers. Simi- 
lar colossal MAE magnitude has been reported for the 
free-standing 4d and 5d TM chains [TT] as well as dimers 
of TMs [ig. Characteristic for the colossal MAE |2] is 
the significant contribution from the variation of mag- 
netic moments and corresponding magnetization energies 
when changing the direction of M, a contribution that is 
absent for conventional 3d magnets. For instance, for 
4x4 W on graphene, this variation of the spin moment 
reaches as much as 0.1/is (Fig. 2a). By looking at the 
MAE values as the function of the band filling and de- 
position density in Fig. 2b, we find that, in contrast to 
the spin moments, the MAE exhbits a much more ir- 
regular behavior as a function of the band filling and 
depends much more sensitively on the adatom density. 
E.g. for W the MAE changes sign depending on the den- 
sity of adatoms, and has a small value in the 4x4 geome- 
try. Such a situation provides a possibility to manipulate 
the magnetization direction with weak external perturba- 
tions, which can be particularly important for transport 
applications as we demonstrate below. 

In search for ability to manipulate the magnetic prop- 
erties of adatoms by external fields, we show in Fig. 2 
the dependence of the spin moments and the MAE of 
W adatoms in 4x4 geometry on the strength of the elec- 
tric field £, applied perpendicularly to the graphene layer 
(along the z-axis). Remarkably, the spin moment of W 
displays a strong dependence on the field strength, espe- 
cially when the magnetization points out of plane. Qual- 
itatively, this response can be characterized by a magne- 
toelectric coefficient a, which relates the change in the 
spin moment to the strength of the f -field in the zero- 
field situation: yLtoA/zs(W) = a£, where is the vacuum 
magnetic permeability constant. For 4x4 W on graphene 
and an out-of-plane magnetization, aj_ amounts to about 
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FIG. 1: Magnetic moments and MAE of 5d adatoms on 
graphene. (a) Magnetic moments due to electron spin, fis, 
calculated without SOC and (b) the MAE of 5d TM adatoms 
on graphene in 2x2, 3x3 and 4x4 superlattice geometry. 
Positive (negative) values of MAE imply an out-of-plane (in- 
plane) easy axis of the spin moments, i.e. perpendicular to 
(in) the graphene plane. 

3 • 10^^'^ G-cm^/V, which is one order of magnitude 
larger than that in Fe thin films [5]. For an in-plane 
magnetization, the magnetoelectric coefficient ay is only 
6 • 10~^^ G-cm^/V, i.e. one order of magnitude smaller 
than a± , which underlines the strong anisotropy of this 
quantity. We observe similarly large variations of the 
magnetoelectric coupling strength in other considered 5d 
TMs. The variation of orbital moment with varying £- 
field is negligible for all systems. 

Striking is the effect of the electric field on the magne- 
tization direction (see Fig. 2b). At zero field {£ = 0) the 
magnetization is in-plane. Applying a negative field of 
magnitude of 0.05—0.4 V/A, values typical in graphene 
field effect transistor structures [55], the sign of the 
MAE is changed. This means that the equilibrium direc- 
tion of the magnetization can be switched from in-plane 
(MAE < 0, for £■ = 0) to out-of-plane (MAE > for 
£ < 0). Supposing that the 5-field is completely screened 
in our system by forming a screening charge Sq^ the vari- 
ation of the MAE with respect to 6q, SMAE/Sq, reaches 
as much as 28 meV/e for negative electric fields, which is 
more than three times larger than that on the surface of 
CoPt slabs [2^, and one order of magnitude larger than 
that of Fe slabs |8| . With a moderate out-of-plane elec- 
tric field of ±0.13 V/A switching of the magnetization 
can be also achieved in 4x4 Hf, and the MAE can be sig- 
nificantly altered in 4x4 Os on graphene (by 10 meV) 
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FIG. 2: Dependence of the magnetic moments (a) and mag- 
netic anisotropy energy MAE (b) of 4 x 4 W on graphene on 
the strength of an external electric field, (a) The spin mo- 
ment of W adatoms fis (in Ms ) as a function of the strength 
of an external electric field £, perpendicular to graphene sur- 
face, for the out-of-plane and in-plane magnetization, respec- 
tively. Negative values of £ correspond to the electric field 
in -|-z-direction, i.e. the direction from graphene towards the 
adatoms. (b) Dependence of the MAE on the strength of an 
external electric field £, sign convention of MAE is consis- 
tent to Fig. 1. Inset displays the difference of spin densities 
Am(r) for the system without electric field and the case with 
a negative electric field of 0.13 V/A. 



and 4x4 Ir on graphene (by « 20 meV). This establishes 
5c? TM adatoms on graphene as a class of magnetic hybrid 
materials with high susceptibility of magnetic properties 
to the electric field, thus making the electrical control of 
magnetism in these systems possible. 

The reason for such a strong magnetoelectric response 
of 5d TMs on graphene can be exemplified for the case of 
W in 4x4 geometry. The local W s- and d-decomposed 
density of electron states without SOC, grouped into 
Ai(s, 42), A3(dj,^,dy2) and A4{d^y,d^2_y2) contribu- 
tions, is presented in Fig. 3a. As we can see, the W 
spin moment of about 1.6/iB originates from two occu- 
pied spin-up and two unoccupied spin-down Ai-states, 
while the exchange-split A3- and A4-states are situated 
above and below Ep, respectively, and almost do not 
contribute. Upon including SOC (e.g. for out-of-plane 
spins) a strong hybridization between the A\, and A3 
and A4 states of both spin occurs around Ep (see also 
Fig. 3b) , which results in a formation of hybrid bands of 
mixed spin and orbital character, while the A| band re- 
mains mainly non-bonding (see Fig. 3c). When an elec- 
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FIG. 3: Electronic structure of 4 x 4 W on graphene. (a) Density of Ai(s,d^2), As{dx 
without SOC. Up and down arrows mark spin-up and spin-down channels, respectively, (b) Band structure without SOC along 
high symmetry lines in the 2D Brillouin zone. Red and blue color of the bands stands for the spin-up and spin-down character, 
respectively. Circles highlight the points at which the gap will open when SOC is considered, (c) Band structures with SOC 
for M||z: without electric field (solid line), with positive (dashed) and negative (dot-dashed) electric field of the magnitude of 
0.13 V/A. The inset shows the anomalous Hall conductivity of the system with respect to the position of the Fermi energy Ef 
for £ = and M||z. (d) Berry curvature distribution of occupied bands in the momentum space (in units of with a as the 
in-plane lattice constant of the 4x4 supercell). The Brillouin zone boundaries are marked with solid lines. 



trie field is applied along the z-axis, the states which 
experience most influence of the corresponding poten- 
tial change are the A|-states, directed perpendicularly 
to the graphene plane. Corresponding modification of 
the band structm'es due to the electric fields can be seen 
in Fig. 3c, where the A| band is shifted downwards (up- 
wards) by negative (positive) applied £ fields, while the 
hybrid bands of mixed character remain almost unaf- 
fected. This causes the redistribution of the electrons 
in the A|-states, its hybridization with the hybrid band 
below the Fermi level and hence the variation of the mag- 
netic moment. This can be visualized by the plot of the 
spin density difference Am(r) for the case with £ = and 
£ = -0.13 V/A (inset in Fig. 2b). It is obvious that, upon 
applying an electric field, a certain amount of spin den- 
sity is transferred from the Ai state of (1^2 character to 
the A3 states. Owing to the difference in the hybridiza- 



tion with the graphene states, seen from the width of the 
corresponding peaks in the density of states, the Ai and 
A3 states have different localization inside the W atoms 
thus leading to a change in the spin moment. 

At last we turn to the prediction of a stable QAHE, 
a manifestation of the quantization of the transverse 
anomalous Hall conductivity. It is motivated by the 
observation that for 4x4 W on graphene with an out- 
of-plane magnetization, A3 and A4 bands of opposite 
spin cross (see circles in Fig. 3b) and hybridize under the 
presence of SOC, forming a global band gap across the 
Brillouin zone (BZ). Thus, 4x4 W becomes an insulator 
upon a spin-orbit driven metal-insulator transition. We 
compute the anomalous Hall conductivity of this system, 
given by a^y = (e^//i)C, with C as the Chern number of 
all occupied bands that can be obtained as a fc-space in- 
tegral C = 5^ f^(k) d'^k. The integrand, fi(k), is the 
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so-called Berry curvature of all states below the Fermi 
level: 

n{k) = X! X! 2Im '^^"'^ ^ ^ ^""^^ ^^""^ ^ '"^ ^ (1) 

where 'i/'nk is the spinor Bloch wavefunction of band n 
with corresponding eigen energy e„k, Vi is the «'th Carte- 
sian component of the velocity operator. The Berry cur- 
vature in reciprocal space, presented in Fig. 3d, displays 
a comparatively complex pattern, with positive contri- 
butions concentrated along the BZ boundary, and large 
negative dips at certain isolated fc-points away from the 
high-symmetry points. Such non-trivial distribution of 
the Berry curvature, which plays a role of an effective 
magnetic field in the fc-space, is typical for complex 
transition-metal compounds |3j. 

The calculated anomalous Hall conductivity as a func- 
tion of the electron filling expressed as energy E relative 
to the Fermi level is presented in the inset of Fig. 3c for 
an out-of-plane magnetization, M||z, and zero f-field. 
We indeed find that the Chern number of all occupied 
states acquires an integer value of +2. According to the 
physics of the quantum Hall effect, this will result in two 
dissipationless and topologically protected edge states on 
each side of a finite graphene ribbon with W adatoms. 
Remarkably, another 32 meV large QAHE gap with the 
C = —2 can be observed at the energy of —0.27 eV be- 
low Ep, originating from the SOC-mediated hybridiza- 
tion between the A3 and A4 bands of the same spin-up 
character, see Figs. 3b and 3d. Such a peculiar situa- 
tion suggests that the topological properties of the QAHE 
state in 4x4 W on graphene can be controlled by tuning 
the position of Ep e.g. via deposition on an appropriate 
substrate. 

More importantly, at equilibrium the magnetization 
direction of 4x4 W lies in-plane, rendering axy = ow- 
ing to the antisymmetric nature of the anomalous Hall 
conductivity with respect to the magnetization direction. 
From our previous observations it follows, however, that 
the direction of M can be conveniently switched to out- 
of-plane by applying a moderate electric field without any 
significant modification of the the electronic structure. 
The gaps remain intact (cf Fig. 3c) , its size unchanged, 
and the value of the anomalous Hall conductivity remains 
quantized at both energies. 

Such non-trivial topological QAHE states occur also 
for all other 5d adatoms (except Pt) on graphene in 
4x4 geometry, and several of them in 2x2 geometry, 
with the QAHE gaps of comparable magnitude, although 
positioned away from the Fermi energy. For example, 
graphene with Re deposited in 4x4 geometry exhibits 
three C = +2 QAHE gaps at -1, -0.15 and -^0.67 eV, 
with the corresponding gap size of 38, 63 and 98 meV, 
respectively. Also, 4x4 Ir and Os on graphene show two 
QAHE gaps (C = +2 for Ir and C = ±2 for Os) at 0.24 



and 0.47 eV for Ir, and -0.38 and 0.15 eV for Os, with 
corresponding gap size in between 11 and 80 meV. The 
fact that several QAHE gaps with different Chern num- 
bers at different energies can occur in the same material 
makes these systems a rich playground for topological 
transport studies. Characteristic of such systems is the 
leading role of the 5d orbitals and strong SOC of heavy 
adatoms in determining the topological properties: as we 
can see from Fig. 3b, the creation of the QAHE gaps can 
occur at the points in the Brillouin zone away from high- 
symmetry points due to hybridization of bands with a 
strong TM character, which are significantly altered by 
the TMs' SOC. This is in contrast to the cases consid- 
ered previously [3 [6l [20] , in which the role of the mag- 
netic adatoms is rather to insert a Rashba spin-orbit or 
exchange fields on the Dirac states at the Fermi energy 
and high-symmetry points in the Brillouin zone. Most 
importantly, however, the strong SOC of 5d TMs leads 
to an increase of the QAHE gap by an order of magni- 
tude, when compared to all QAHE systems suggested be- 
fore, guaranteeing a strong topological protection against 
defects, structural disorder or thermal fluctuations, and 
thus opening a route to QAHE at room temperature. 

To summarize, we predict that graphene decorated 
with 5d transition-metals is a hybrid material display- 
ing remarkable magnetic and transport properties, which 
may be utilized experimentally for the purposes of topo- 
logical transport applications at reasonable temperatures 
and conditions. The fact that heavy transition-metals de- 
posited on graphene display very stable magnetism due 
to colossal magneto-crystalline anisotropics with exciting 
possibilities in engineering a robust quantum anomalous 
Hall effect with a large band gap stems from the strong 
spin-orbit interaction of these elements. The extreme 
sensitivity of these hybrid systems to an electric field that 
we predicted is likely to be enhanced when graphene is 
deposited on a ferroelectric substrate. Thus 5(i-graphene 
hybrid materials provide a unique functionality in tai- 
loring the required magnetic and topological properties, 
not likely to be achieved in conventional magnetic ma- 
terials based on 3d transition-metals. Considering the 
breadth of this material class including other possible 
heavy adatoms, such as the Ad TMs or 4/-elements, we 
anticipate that room temperature QAHE may become 
possible in these systems and we encourage intense ex- 
perimental research on these materials. 

Method 

The theoretical investigations are based on density 
functional theory |28j . which is 'first-principles' in that 
it requires no experimental input other than the nu- 
clear charges. We apply the generalized gradient ap- 
proximation |29j to the exchange-correlation potential 
and use the full-potential linearized augmented plane 
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wave method (FLAPW) as implemented in the FLEUR 
code [30]. The self-consistent calculations with SOC were 
carried out with the cut-off parameter k^ax of 4.0 bohr~^ 
and 144 /c-points (for the 4x4 cases) in the full two- 
dimensional Brillouin zone. The muffin-tin radius of 
2.7 bohr (1.2 bohr) was used for TM (carbon) atoms. 
To calculate the transverse anomalous Hall conductivity 
accurately, we applied the Wannier interpolation tech- 
nique, using the FLEUR and WannierQO codes to construct 
the maximally- localized Wannier functions [5TH55] . For 
all systems, the distance between the TM atoms and the 
rigid graphene layer was structurally optimized by total 
energy minimization neglecting SOC. We confirmed that 
the relaxations of the carbon atoms within the graphene 
layer are indeed very small, and neglecting them does not 
lead to any changes in the calculated quantities. TMs 
bond covalently with graphene [2T] , thus we suspect that 
the van der Waals interaction plays a minor role. 
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